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Abstract 

Semiclassical (stochastic) wave equations are proposed for the cou- 
pled dynamics of atòmic quantum states and semiclassical radiation 
field. All relevant predictions of Standard unitary quantum dynamics 
are exactly reproducible in the framework of stochastic wave equa- 
tion model. We stress in such a way that the concept of stochastic 
wave equations is not to be restricted to the widely used Markovian 
appr oximation . 
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1) Introduction. Since many years it has been understood that one can 
not reproduce all statistical predictions of quantum theory from pure clas- 
sical statistics One can do it for certain predictions after all. Reduced 
dynamics of quantum subsystems (i.e. effective dynamics of open quantum 
systems) can always be interpreted as hybrids of quantum and classical mech- 
anisms. The corresponding semiclassical dynamics (which may on equal 
footing be called semiquantized) is described by stochastic wave equations 
(SWEs). Since ten years they have been attracting permanent interest in 
research of foundations and, since the 90's ||, in quantum òptics, too. 

Consider a system of àtoms interacting with radiation field. Basically, 
all attainable results are approximations: SWEs have been derived in weak 
coupling Markovian limit Q|. In the present work, however, we investigate 
the exact unitary dynamics of the atom+radiation system and we ask the 
following question. Can we describe the radiation field in terms of classical 
(stochastic) variables instead of quantum ones in such a manner that all rel- 
evant predictions remain exactly identical with those of the unitary quantum 
theory? The answer will be affirmative in almost all respects. Forthcoming 
results were anticipated once f| by detailed relativistic calculations. Here we 
intend to give a short account on exact SWEs. Readers interested in practi- 
cal rather than conceptual aspects may, after reading Section 2, undertand 
SWEs directly from Section 5. 

2) The quantum model. For simplicity's sake, we consider single mode 
(cavity) quantum-electrodynamics. The Hamiltonian for the composite sys- 
tem of the atom plus radiation field has the following general form 

H = H ato + u~a)à + Ta) + Pa (1) 
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where H ato is the atòmic Hamiltonian, à is the absorption operator for the 
single radiation mode of frequency u. In the interaction terms, J is atòmic 
operator proportional to charges and currents. As usual, we assume that at 
time t = the quantum state of the system is a product of the atòmic and 
cavity radiation states, respectively: 

*o = W to V ad - (2) 
In interaction picture, the quantum state \l/ evolves unitarily: 

-iÈ!(t)%, (3) 



dt 

with interaction Hamiltonian 

È I (t)=j(t)iï + H.C. (4) 

where j(t) = e lu}t J(t) is the "rotated" version of the interaction picture cur- 
rent J(t). Though the present work concerns pure state equations it will 
nonetheless be convenient to define a density operator p = too. 

3) Stochastic field variables, semiclassical states and observables. The bà- 
sic field variables à, a' do not commute: [a, a'] = 1. We shall, nevertheless, 
establish a certain natural correspondence between them and their classi- 
cal counterparts a,a*. Let us introduce a special notation for symmetric 
products: 

à c Ò = -{àÓ + Òà) = {à, 0} (5) 

where Ò is arbitrary operator and à c is called " superoperator" (c.f. p and 
references therein). The superoperators à c , commute with each other 0: 

[a c , à\]Ò = {a, {at, Ó}} - {at, {a, Ó}} = 0, (6) 



hence we shall make the natural correspondence, mentioned above, between 
the classical complex variable a (or a*) and the superoperator à c (or aj). 

We describe the radiation field in terms of the classical stochastic variables 
a, a* coupled to the quantized atòmic ones. The state of the quantum system 
will then be represented by the semiclassical "density": 

p(a, a*) = tr rad (ó(a - à c )5(a* - à\)pj , (7) 

which is expected to be density operator for the atom and phase space distri- 
bution for the field simultaneously. It should satisfy the positivity condition 
p(a, a*) > for all complex vàlues of a, a*. Then tr ato p(a, a*) is the reduced 
phase space distribution p rad (a, a*) of the radiation field mode and, alter- 
natively, / p(a, a*)da*da is the reduced density operator p ato of the atom. 
Similarly to Eq. (7), we define semiclassical counterpart of Hermitian ob- 
servable F by: 

F(a, a*) = tr rad (5(a - a c )5(a* - àl)P) . (8) 

We do however not consider all Hermitian observables but those where the 
correspondence is invertible: 

tr(F{à c ,al)Ò) =tr(FÓ) (9) 

for all Ò. This includes a pretty large class of Hermitian operators. 

The equivalence of the semiclassical picture with the full quantized one 
relies upon the following equation: 

tr (P(à cj àl)pj = tr ato j F(a, a*)p(a, a*)da*da (10) 

which follows from Eqs. (7-9) and from the commutativity of à c and aj. 

3 



4) The positívity issue and its solution by coarse-graining. Let us write 
the phase space distribution of the radiation field, defined by tracing Eq. (7) 
over the atòmic states as well, in Fourier-representation: 

p - (a , a . ) = J_/exp(-Aa· + A·a) í r(exp(-A·a c + Aa:)p)<ÍA· íí A. (11) 

By substituting the identity 

exp(— X*à c + Xàl)p = exp(— A*a/2) exp(Àày2)/3exp(Aa''/2) exp(— A*à/2) 

(12) 

which follows from the rules (5,6), one can apply the identity 
exp(i) exp(B) = exp([i, B}/2) exp(i + B) 

to obtain 

1 I" 

p rad (a } a*) = — J exp(-Aa* + X*a)tr (exp(-A*a + Aa f )p) d\*d\. (13) 

One recognizes that p rad (a, a*) is identical to the well-known Wigner-function 

introduced long ago 0. However, the Wigner-function can not be inter- 

preted as phase space density since it may take negative vàlues as well. From 

Husimi's work we know that a minimum coarse-graining will cure the prob- 

lem ||. Consequently, we apply a minimum coarse-graining to the "sharp" 

distribution p(a,a*) defined by Eq. (7). Let a ,a^ be noisy complex field 

variables of vacuum distribution (2/7r) exp(— 2|ao| 2 ). Coarse-grained semi- 

classical state thus reads: 

- r exp(-2|a | 2 ) _ * *v, *, , nA s 

p{a,a) = / — p(a + ao,a + a )da dao. (14) 

J 7T/2 

The corresponding coarse-grained phase space density p rad (a, a*) = trp(a, a*) 
is identical to the so-called Q-function of the radiation field (see, e.g., in M). 



As is well-known, the Q-function is l/7r-times the diagonal element of the 
density operator between coherent states. Accordingly, the corresponding 
relation holds true for p(a, a*) itself: 

p(a, a*) = — (a, a*\p\a, a*) (15) 

71 

where \a,a*) is coherent state of the radiation field mode. [Two remarks 
on our notations are in order. First, we indicate the explicit dependence 
of coherent states on both a and a* . Second, the quadratic form above has 
to be understood on the factor Hilbert-space of the cavity while the result- 
ing expression is still an (a, a*)-dependent density operator in the atom's 
Hilbert-space.] The form (15) guarantees the positivity of the coarse-grained 
semiclassical density. 

One also introduces coarse-grained semiclassical observables: 

F(a,a*) = í — ^j a °l - p( a + a , a* + a* )da* daQ. (16) 

J 71/2 

Then a coarse-grained version of the "sharp" equivalence equation (10) fol- 
lows: 

tr ^F(a c , à[)p^ = tr ato j F(a, a*)p(a, a*)da*da, (17) 

i.e. the statistics of coarse-grained quantum observables are still exactly 
reproducible from the coarse-grained semiclassical state p(a,a*). 

One would worry because coarse-grained distributions are believed to 
loose information. In typical quantum òptics applications things prové not so 
bad as long as one considers expectation vàlues of polinoms of the absorption 
and emission operators. Then the additional noise is easy to subtract. Let us 
assume that we are interested in the expectation value of the photon number 
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operator à'à in the current quantum state p and we want to calculate it from 
p(a,a*). It is plausible to choose F(a,a*) = \a\ 2 . Then, via Eq. (16), let us 
calculate the corresponding coarse-grained observable F(à c , aj). We obtain 
àtà c = à\à c + \. The equivalence condition (17) thus leads to: 

tr ({àlà c + ^)pj = J \a\ 2 p rad (a, a*)da*da. (18) 

Due to the identity tr {{a[à c + \)pj = tr{à^àp) + 1, Eq. (18) yields: 

tr(iïàp) = í \a\ 2 p rad (a,a*)da*da-l. (19) 

This is the simplest example to illustrate how exact quantum expectation 
vàlues are reproduced from the coarse-grained semiclassical state p(a,a*). 

5) Exact stochastic wave equation. Remind that the composite system of 
the atom+radiation was originally assumed to be in a pure quantum state 
* '. Let us introduce a coherent state representat ion of 

V>(a,a*) = -^(a,a*|*>. (20) 
V 71 " 

The scalar product on the RHS is to be taken on the factor Hilbert-space of 
the cavity. [The author has failed to find a more suitable compact notation, 
see also the second remark after the Eq. (15).] Strictly speaking, ip(a, a*) is 
wave function of the field mode coordinates a, a* while, on the other hand, 
it is state vector (of arbitrary representation) in the atòmic Hilbert-space. 
It follows from Eq. (15) that the wave function (20) just yields the coarse- 
grained semiclassical density in the form 

p(a, a*) = ip(a, a*)tp\a, a*). (21) 

So, it is natural to call ip(a,a*) semiclassical wave function. Let us obtain 
its equation of motion. 
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In coherent state representation the following operator correspondences 
should be used (c.f. JTTJ ) : 

a*ip(a, a*), f - + ^- J ^(a, a*). (22) 

V 2 aa* y 

Hence the Hamiltonian equation of motion (3) leads to the following equation 
of motion for the semiclassical wave function: 

j t = -ï (a j{t) + -f(t) J ip t {a, a ) - if(t) — — . (23) 

This wave equation, combining unitary and stochastic dynamics, is our 
central result. By means of this equation, exact predictions of the unitary 
quantum theory can be otained JTTJ]. 

Let us summarize the SWE method. Assume the atòmic state at t — 
is i/jQ to while the cavity's state is the vacuum state |0, 0). Then, according to 
Eq. (2), the semiclassical (stochastic) wave function (20) initially takes the 
form 

V>o(a,a*)=Vr eXP( Í^ |2/2) - (24) 



l'K 

Let us switch on the interaction. Then the conditional atòmic wave function 
^>t(a, a*) evolves in function of the state (a, a*) of the classical field, according 
to the Eq. (23) while the (coarse-grained) probability distribution of the 
classical field variables a, a* is given by 

p t (a,a*) = ||^ t (o,o*) f. (25) 

The expectation value of an arbitrary semiclassical observable F(a, a*) will 
exactly reproduce the quantum expectation value of the corresponding coarse- 
grained Hermitian observable, i.e. 

tr ^F(à c , à],)ptj = j ' ipj(a,a*)F(a,a*)?p t (a,a*)da*da. (26) 
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where F(a, a*) is derived from F(a, a*) by the minimum coarse-graining (16). 

The existence of exact SWE was pointed out in Ref. || in context of the 
relativistic quantum-electrodynamics. That time, however, the stochastic 
field was not identified as semiclassical radiation field. A generalization of 
our SWE (23) for infinite number of field modes being initially in thermal 
equilibrium states at nonzero temperature seems straightforward. 

6) Outlook. The main result of the present work is a suggestion that 
semiclassical SWEs could exactly reproduce the quantum physics of both 
atòmic and radiation degrees of freedom. In our opinion, the very existence 
of exact stochastic models is a novelty of fundamental interest. Conceptual 
elements of equivalence with unitary quantum theory are symmetrized prod- 
ucts of non-commuting operators and a minimum coarse-garining of states 
and observables. Our SWE establishes the possibility of making eflicient 
Monte-Carlo simulations in the non-Markovian regime since using classical 
variables for the radiation field represents a fair simplification. Of course, 
the re-derivation of the extensively used Markovian SWEs is possible and 
desirable. 

This work was supported by the grants OTKA No. 1822/1991 and 
T016047. 
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